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INFINITE RATE MUTUALLY CATALYTIC BRANCHING^ 

By Achim Klenke and Leonid Mytnik 

University of Mainz and Technion Haifa 

Consider the mutually catalytic branching process with finite 
branching rate 7. We show that as 7 — > cxj, this process converges 
in finite-dimensional distributions (in time) to a certain discontin- 
uous process. We give descriptions of this process in terms of its 
semigroup in terms of the infinitesimal generator and as the solution 
of a martingale problem. We also give a strong construction in terms 
of a planar Brownian motion from which we infer a path property of 
the process. 

This is the first paper in a series or three, wherein we also con- 
struct an interacting version of this process and study its long-time 
behavior. 

1. Introduction and main results. 

1.1. Motivation. In [5], Dawson and Perkins introduced a population 
dynamic model of two populations that live on a countable site space S. 
The individuals migrate between sites and, at any given site, perform a 
critical branching process with a branching rate proportional to the local 
size of the population of the respective other type. 

More precisely, Dawson and Perkins considered the system of coupled 
stochastic differential equations (SDEs) (taking nonnegative values) 

(1.1) dY,^t{k) = {AY,^t){k) dt + ^-fYi^tik)Y2,tik) dWuik), i = l,2,kGS. 

Here, A{k, I) = a{k, I) — Ij^j (/) is the (7-matrix of a Markov chain on S with 
symmetric jump kernel a, {Wi{k),k £ S,i = 1,2) is an independent family of 
Brownian motions and 7 > is a parameter. 

Dawson and Perkins showed that there exists a unique weak solution 
of this SDK taking values in a suitable subspace of ([0,oo)^)'^ with some 
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growth condition. Furthermore, this process is a strong Markov process. 
While existence of a weak solution is rather standard due to the procedure 
proposed by Shiga and Shimizu [16], weak uniqueness was shown using a 
certain self-duality of the process established in [13]. We will describe the 
duality in detail below, in (2.4). 

A main result of Dawson and Perkins is a dichotomy in the long-time 
behavior of the solutions depending on whether A is recurrent or transient 
(assuming some mild regularity condition on A). For recurrent A (fulfilling 
the regularity assumption), the types segregate, while for transient A, there 
is coexistence of types. More precisely, let 

k£S 

denote the total mass processes {i = 1,2) and assume that Mifi,M2fl < oo. 
Then Mi and M2 are continuous orthogonal nonnegative L'^-martingales. 
Let Mj^oo = hmt^oo ^i,t denote the almost sure limit. Dawson and Perkins 
show that E[Mi,ooM2,oo] = if ^ is recurrent and E[Mi,ooM2,oo] = Mi,oM2,o 
if A is transient. Furthermore, in the recurrent case, the joint distribution of 
(Afi,oo; Af2,oo) equals Q(Mio,a^2o)' where, for x G [0,oo)^, Qx is the harmonic 
measure of planar Brownian motion in [0, 00)^. That is, if B = (i?i,i?2) is a 
Brownian motion in started at x and r = mi{t > 0: Bt ^ (0, 00)^}, then 
Qx is the probability measure on 

E:= [0,oo)2\ (0,oo)2 

given by 

(1.2) Qx = Px[Br£-]. 

The explicit form of the densities of Qx can be found in (2.5). 

Via the self-duality of the mutually catalytic branching process, its total 
mass behavior for finite initial conditions provides information on the local 
behavior if the initial condition is infinite and sufficiently homogeneous. For 
X G [0,00)^, let X denote the state in ([0,00)^)*^ with x^{k) = Xi for all fc G 5, 
i = l,2. Assume that Yq = x. Then 

lim P,[Fi,t(0)F2,t(0) >0] >0, 

" A is transient, that is, types can coexist locally. On the other hand, for 
recurrent A, the distribution of Yt converges weakly to J SyQx{dy), that 
is, to a spatially homogeneous point y, where y is sampled according to the 
distribution Qx- Hence, in the recurrent case, the two types segregate locally 
and form clusters. The assumption that the initial point is constant can be 
weakened to an ergodic random initial condition (see [3]). 
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The starting point for this work was the wish to obtain a quantitative 
description of the cluster growth in the recurrent case. We wiU only briefly 
describe the heuristics. Dawson and Perkins also constructed a version of 
their process in continuous space M instead of S as the solution of a stochastic 
partial differential equation 

dY t(r) I • 

(1.3) -^j^ = ^Yi^t{r) + ^lYr^t{r)Y2,t{r)Wi{t,r), rG]R,i = l,2, 

where Wi and W2 are independent space-time white noises and A is the 
Laplace operator. As A on M is recurrent, types also segregate here. Now, 
due to Brownian scaling, if we denote by Y^ the solution of (1.3) with that 
given value of 7, then we obtain 

(1.4) V^[{Y^{rVT)\^^ G •] = V^[{Yf{r)U^ E •]. 

Equation (1.4) shows that clusters of Y\^t grow like y/T and that a better 
understanding of the precise cluster formation can be obtained by letting 
7 — )• 00 for fixed time. Hence, we aim to construct a model X, that is, in 
some sense, the limit of Y'^ as 7—7-00. 

In this paper, we construct X in the simple case where S" is a singleton 
and where the migration between colonies is replaced by an interaction with 
a time- invariant mean field. This is a first step toward the investigation of 
the model involving infinitely many sites. We give characterizations of the 
process X via an infinitesimal generator, as the solution of a well-posed mar- 
tingale problem and as the limit of Y^ as 7—7-00. Finally, we give a strong 
construction of the process via a time-changed planar Brownian motion. 
This will also serve to derive path properties. 

In two forthcoming papers, we construct the infinite rate process on a 
countable site space S via a stochastic differential equation with jump-type 
noise and give a characterization via a martingale problem [9]. Furthermore, 
we will investigate the long-time behaviour and give conditions for segre- 
gation and for coexistence of types [10]. An alternative construction via a 
Trotter product approach is carried out in [11] and [14]. 

1.2. Results. We now describe the one-colony process which is the sub- 
ject of investigation of this paper. Assume that S is a singleton and that im- 
migration and emigration come from and go to some colony that is thought 
to be infinitely big and whose effective population size (for immigration) 
is € [0,00)^. Furthermore, let c > be the rate of migration. Hence, we 
consider the solution Y = Y"*''^'^ of the stochastic differential equation 

(1.5) dY,,t = c{e^ - y,,t) dt + yJlYi^tY2,t dWi,t, i = l,2. 

This model can be thought of as a version of the model defined in (1.1) where 
the migration between colonies is replaced by an interaction with a time- 
invariant mean field 9 or with an infinitely large reservoir whose types have 
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proportions 9i and ^2- (In fact, in [2] it was shown (Proposition 1.1) that 
y7,c,6» ^j-jggg McKean-Vlasov hmit of solutions of (1.1) with symmetric 

interaction on a complete graph S.) More formally, the interaction term AY 
is replaced by a drift c{9i — 5^,t). It is this simplification of the interaction 
that allows for a tractable exposition in this article. Note that as t — )• oo, the 
process without drift (c = 0) converges almost surely to some random x £ E. 
Hence, in the case c = 0, if we let 7 — t- 00, then the limiting process would 
be trivial: if it starts at x G E, then it stays at x forever. See Section 2 for a 
more detailed description of the process Y solving (1.5) (finite 7 process). 

On a heuristic level, as the stochastic term in (1.5) defines an isotropic 
two-dimensional diffusion, that is, a time-transformed planar Brownian mo- 
tion, if we let 7 — )• 00, then we should end up with a process where the 
stochastic part is a planar Brownian motion at infinite speed, stopped when 
it reaches the boundary of the upper-right quadrant. That is, the limiting 
process X should be a Markov process with values in E. When x is the cur- 
rent state and the drift moves it to x -|- c{6 — x) dt, this point should instanta- 
neously be replaced by a random point chosen according to Qx+c{e-x)dt- 

We 

will, in fact, be able to describe this infinitesimal dynamics both in terms 
of a martingale problem and in terms of a generator of Markov transition 
kernels. However, we first define X via an explicit transition semigroup and 
show that it is the limit of Y'^'''^'^ as 7—7-00. Let 

(1.6) Ci{E) :=\f:E-^Cis cont. and lim f{u,0) = lim f{0,v) is finite) 

L )-00 >00 J 

equipped with the supremum norm ||/||oo = sup^.g£;|/(x)|. 

Definition 1.1. Let c > and 9 G [0, 00)^. For t>0 and x £ E, define 
the stochastic kernel pt by 

Pt{x, •) ■=p1'^{x, •) := Qe-ct^+(i_e-ct)0. 

Define the contraction semigroup S = {St)t>o on Ci{E) by 

Stf{x)= I f{y)ptix,dy). 
Je 

The Markov process X = X'^'^ with state space E, cadlag paths and tran- 
sition kernels {pt)t>o is called the infinite rate mutually catalytic branching 
process (IMUB) with parameters {c,9). 

In order for this definition to make sense, in Proposition 3.2, we will show 
that {St)t>o is, in fact, a Markov semigroup. 

Proposition 1.2. X'^'^ is a Feller process and has the strong Markov 
property. It is ergodic and the unique invariant measure is Qg. 
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Proof. The map x i— )■ Qx is continuous, hence x i— ?• P( (x, •) is also contin- 
uous, that is, X'^'^ is a Feher process. Since Qx = 5x for x ^ E, the semigroup 
S is strongly continuous. Hence, by the general theory of Markov processes, 
there exists a cadlag version of X that is strong Markov (see, e.g., [15], 
Chapters III. 7 and 8). 

Ergodicity and the explicit form of the invariant measure are trivial. □ 

Theorem 1.3 (X"^'^ as an infinite rate process). Assume that Y^''^'^ = 

c 

Xq = x £ E for all 7 > 0. As 7 — ?• 00, the finite- dimensional distributions 
ofY'^''^'^ converge to those of X'^'^ . 

Note that in Theorem 1.3, trivially, we do not have convergence in the 
Skorohod path space, since continuous processes do not converge to discon- 
tinuous processes in that topology. 

In addition to the convergence of the finite-dimensional distributions, we 
also have convergence of the pth moments for p G [1; 2) [but not for p = 2, of 
course, since for x £ (0,oo)^, the measure Qx does not possess finite second 
moments, as can be easily derived from its density formula (2.5)]. Hence, on 
a suitable probability space, we have L^-convergence of y^''^'^ to X'^'^ . 

Theorem 1.4 (L^'-convergence). Assume that 5^''^'^ = X^^ = x e E for 
a//7>0 and letpG [1,2), t>0. 

(i) For every 7 > and i = l,2, we have 

E.[(y,7'^)n<E.[(X^;f)^]<oo. 

(ii) On a suitable probability space, for i = 1,2, we have 

It can be seen from the proofs of Theorems 1.3 and 1.4 that the statements 
of these theorems also hold for Yq''^'^ = x £ [0,oo)^ and t > if we replace 
Xq^ by a random point chosen according to Qx- 

Remark 1.5 (Trotter product approach). While in the one-colony case 
considered in this paper, it is easy to explicitly write down the semigroup for 
the infinite rate mutually catalytic branching process X'^'^ , it is less obvious 
how to construct an interacting version of the process on a countable site 
space. One possibility is the Trotter product approach that is used in [11] 
and [14]. Here, we briefly sketch it for X'^'^ . 

In the classical setting, the Trotter product approach works as follows. In 
order to construct a solution y^'^^'^ of (1.5), in time intervals of length e, one 
could alternate between a solution of the pure drift equation (7 = 0) and the 
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pure stochastic noise equation (c = 0). As e 4, 0, this process converges to a 
solution of (1.5). 

If we let 7 —7- 00, then the noise term results in an instantaneous jump to 
a point in E chosen according to Qy, where y is the value of Y at the end of 
the preceding "drift interval." More formally, let {^{k,x),k G N,x G [0,oo)^) 
be an independent family of -E-valued random variables with distribution 
x)] = Qx - For t S [ke, {k + l)e), let be the solution of the differential 
equation 

dXf = c{e - Xt) dt, 

that is, 

Xf = e"<*-''^'^ Xl^ + (1 - e-^^^-^'^)e. 

Let 

^ffc+iv_:= lim Xt = e-'^Xf^^ + {l-e-^)e 

and define 

^lk+l)e = '^(^ + li^ffc+l)e-)- 

One can prove that X^ converges in distribution in the Skorohod topology 
on the space of cadlag paths to X'^'^ (see [11] and [14]). 

While, in Definition 1.1, we gave an explicit formula for the transition 
kernels of X, it is also interesting to characterize the process X via its in- 
finitesimal dynamics. In Section 5, we investigate the generator Q of the 
semigroup S. For a certain class Cf{E) C Ci{E) of smooth functions / (see 
Definition 5.1), we give an explicit formula for ^/ as an integro-differential 
operator. Using the classical Hille-Yoshida theorem, we show that the re- 
stricted operator Q = Q\(j2(^^-^ uniquely defines {St)t>o (Theorem 5.3). Fur- 
thermore, we show that Q restricted to an even smaller space V of functions 
that appear in the duality for X still uniquely defines the process X via a 
martingale problem (Theorem 5.4). To define it is crucial to study (for 
suitable functions /) the limit 

lim*~n'5t/(2;) - f{x)) = Irni e^^ /dQ^+ede-x) - f{x)^ , 

which will also clarify the jump structure of the process X. The description of 
the exact form of the operator Q and the precise statements of the theorems 
are a bit technical, so these are deferred to Section 5. 

While, for Proposition 1.2, we used general construction principles of 
Markov processes, here, we provide an explicit strong construction of the 



INFINITE RATE MUTUALLY CATALYTIC BRANCHING 



7 



process X in terms of a given planar Brownian motion B. This construction 
also allows certain path properties to be investigated. 
Assume Bq = Q. For z S M^, we write 

[2:,00) = [Zl,00) X [2:2,00) 

for the rectangular cone northeast of z. For x G [0,oo)^, let 

(1.7) T^:=\ni{t>Q:Bti[-x,oo)} 
and 

(1.8) D^:=Br,+xeE. 

For x,y £ M?, we write y < x if x G [y,oo), that is, if yi < xi and 1/2 < X2- 
For x G [0,00)^, we define the a-algebra 

(1.9) :F^ = a{Dy:y<x). 

In Lemma 3.1, we will show that D is a Markov process with respect to 
(•^i')a;G_[0,oo)2- 

Let 9: [0,oo) — [0,oo)^ and c: [0,oo) — [0,oo) be measurable and locally 
integrable. For < s <t, define 



(1.10) C{s,t) = exp(^- c{r)dr^ and E{s,t) = J^' 



C(0,r 



dr. 



Theorem 1.6. Let x e E and define the process X'^'^ by 

= C(0,t)Z?,+s(o,t), t>0. 

Then X'^'^ is a time-inhomogeneous Markov process on E with cddldg paths 
and with transition probabilities 

(1.11) Ps,t{^r)=Qcis,t)z+C{0,t)E{s,t) forO<s<t,zeE. 
In particular, for 9 = 6 £ [0, 00)^ and c = c>0, 

(1.12) Xf'''^ = e"''*L>^+(ect_i)0 

is an infinite rate mutually catalytic branching process with parameter {c,6), 
see Figure 1. 

It is tempting to use this strong construction of X'^'^ in order to define an 
interacting version of the infinite rate mutually catalytic branching process 
on a countable site space S, where c6k{t) at site k £ S reflects the migration 
from neighboring sites to k. However, in this paper, we do not pursue this 
topic. Rather, we use the strong construction in order to derive a path 
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property of X'^'^ via a result of Le Gall and Meyre [12] on the cone points 
of planar Brownian motion. 

Recall that a measurable set Ac E is called polar for X'^'^ if for all x & E, 
we have 

Px[Xt'^ G A for some t>0]=0. 



Theorem 1.7. The point OgE is polar for X'^". 

1.3. Organization of the paper. In Section 2, we give a detailed descrip- 
tion of the duality for the process with finite branching rate. In Section 
3, we establish a similar duality for the infinite rate process and use it in 
order to show the convergence in Theorems 1.3 and 1.4. In Section 4, we 
justify the strong construction of Theorem 1.6 and also prove Theorem 1.7. 
Finally, in Section 5, we describe the infinite rate process in terms of its 
infinitesimal dynamics and state and prove the theorem on the construction 
via the Hille-Yoshida theory (Theorem 5.3) and via a martingale problem 
(Theorem 5.4). 




Fig. 1. Strong construction of X^^'^'^^'^^ with Xo = x = (0,1) via a planar Brownian 
motion. Here X^^'^^'^' = e-*''^((0, 1) + &t + (2, l)(e*''2 - 1)) for t = 0,l, 2, 3. 
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2. Duality of the finite 7 process. A major tool for the investigation of 
mutually catalytic branching processes is a self-duality for the process. As 
it turns out to be crucial also for the limiting case of infinite branching rate 
(7 = 00), we describe this duality here in more detail. For x = (xi,X2) and 
y = (2/1)2/2) G K^, we introduce the lozenge product 

(2.1) xoy:= -{xi+X2)iyi+y2) - X2)iyi - ^2) 
(with i = and define 

(2.2) F{x,y) = exp{xoy). 

Note that xo y = y o x. Furthermore, define the "scalar product" 

(2.3) {x,y) =xiyi+X2y2 for x, y G [0, 00)^. 
For X = {x{k))k&s and y = {y{k))kes, we write 

H{x, y) = exp ( ^ x{k) o y{k) 

If Y is the process defined in (1.1) started in state y and Y is the process 
started in some suitable y (such that all sums are finite), then the duality 
reads (see [13], equation (2.5)) 

(2.4) ^y[H{YtM='E^y[H{yM- 

In fact, this duality also holds for asymmetric ^ if y is a solution of (1.1) 
with A replaced by its transpose A* . As this mixed Laplace and Fourier 
transform H is measure determining ([13], Lemma 2.5), the duality yields 
uniqueness of the solutions of (1.1). Furthermore, it provides a tool for trans- 
lating local properties of the solutions into global properties and vice versa. 
If X = (u,f) £ (0,00)^, then the harmonic measure Qx [recall (1.2)] has a 
one-dimensional Lebesgue density on 

E ■= ([0, 00) X {0}) U ({0} X [0, 00)) 

that can be computed explicitly 

4 uvu -f-n 

■ au, if f = U, 



4 uvv -f-n 

■ aw, if u = U. 



(2.5) Q{u,v){d{u,v)) = 

. TT Au^v"^ + (u^ — v'^Y 

Furthermore, trivially we have 

(2.6) Qx = 5x if a; G E. 

We now turn to the situation of only one colony. We consider the solution 

Z = (Zi,Z2) of 

(2.7) dZ,,t = ^/lZi^tZ2,tdWi,u i = l,2, Zo = zG[0,oof. 
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By Theorem 1 of [4], there is the unique strong solution to the above equa- 
tion. 

Clearly, Z\ and Z2 are orthogonal L^-martingales and hence they con- 
verge almost surely to some random variable Z^o = (-Z^i,oo) ■^2,00)- As Z is 
an isotropic diffusion on [0,oo)^, it is a time-transformed Brownian motion. 
Thus Zoo has the same distribution as a planar Brownian motion B started 
at z and stopped (at time r) upon leaving (0,oo)^, that is [see (2.5)], 

l^z\Zcx\ = ^z[Bt] = Qz- 

(We denote by Cx[Xt] = Px[Xt e •] = P[Xt £ • \ Xq = x] the distribution of 
the process X at time t when started at x.) It is easy to see that, in fact, 

■.= mf{t>0:Zt G E} < 00 almost surely, 

and that 

Zt = Z^z for all t>T^. 

Clearly, increasing 7 amounts to speeding up the process. Hence, in the 
limit, we would have a process that instantaneously jumps from z to a 
random point (picked according to Qz) and then stays there. In order to ob- 
tain a more interesting limiting process, and with a view toward interacting 
colonies, we introduce a drift term and consider the equation (which was 
analyzed in more detail in [2], Propositions 1.1 and 1.2) 

(2.8) dY.^t = c{e^ - Yi^t) dt + V7n,t^2,t dWi^t, i = l,2. 

Here, c > and 9 S [0,oo)^ are parameters of the process. It is standard to 
show that (2.8) has a weak solution. Weak uniqueness can be obtained via 
duality. We first outline the general picture for the duality that comes from 
the interacting colonies case and then give an explicit computation for our 
special situation. 

Let us consider a two-colonies model with site space 5 = {1, 2}, where Y 
is the size of the population at site 1 and the size of the population at site 
2 is constant and equals 9. This amounts to a migration matrix 

(2.9) .4=(-; 

and to branching rates 7(1) = 7 (at site 1) and 7(2) = (at site 2). Note that 
the approach of Dawson and Perkins does not require that the branching rate 
be constant; neither does it require that the migration matrix be symmetric 
or a (7-matrix. (At least if S is finite — otherwise, certain regularity conditions 
have to be imposed.) Dawson and Perkins use a duality with respect to a 
process Y with migration matrix A* (the transpose of A) and with the same 
branching rates as Y to show weak uniqueness of Y . 
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Let us now construct the dual process explicitly. We will later use this 
approach in order to construct a dual for the 7 = oo limiting process. Let 
y = (y(l),y(2)) € ([0,oo)^)^ and let Z be the unique strong (by Theorem 1 
of [4]) [0, 00) ^-valued solution of 

(2.10) dZ.^ = ^jZi^tZ2,tdW^,t, i = l,2, Zo = y{l). 
Define a process Y on ([0,oo)^)^ by 

(2.11) Yt{l) = e-''Zt and Yt{2) = y{2) + [ ce-^"- Zr dr. 

Jo 

Note that this y is a solution of (1.1) with 5 = {1, 2}, with site-dependent 
branching rate 7(1) = 7, 7(2) = and with A from (2.9) replaced by A*. 
In particular, y is a time-homogeneous Markov process. We also get the 
time-homogeneous Markov property via an explicit computation: 

Yt+s = [e-''^'^'^Zt+s,y{2) + fj' ce-^'Zrdr^ 

= (^e-''{e-''Zt+s),m) + 1^ ce-'''Zrdr + l\e-'''{e-''Zt+r)dr^ 

= (e— Z;,y'(2)+^'ce""Xrf^), 

where Z'^ = C'^Zt+r and y'{2) = Yt{2) = y{2) + ce'"'' Zr dr. Clearly, Z' has 
the distribution of a solution of (2.7) with y'(l) := Z'q = Yt{l). 
For x,x',y,y' G [0,oo)^, recall that 

(2.12) H{{x, x'), {y, y')) = F{x, y)F{x', y'). 



Proposition 2.1 (Duality). LetY andY be defined by (2.8) and (2.11), 
respectively. Then, for all y G [0, 00)^, y G ([0,oo)^)^ and t > 0, we have 

(2.13) ^y[H{{Yt,e),y)]=Y:.y[H{{y,e),Yt)]. 

In particular, if Z is a solution of (2.10) with Zq = z ^ [0,oo)^, then 



(2.14) EjF(yt,z)]=E, 



F{y,e-''^Zt)F(e, ce-'^^Z^dr 



A similar duality was derived in [2], Lemma 4.2. Before we prove the 
proposition, we have to collect some properties of the derivatives of F. We 
omit the proof of the following lemma. 
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Lemma 2.2 (Derivatives of the duality function). Denote the partial 
derivatives of F by 

ViF(x, y) := ^F{x, y), V2F{x, y) := y) 

and define the Laplace operators Ai and A2 hy 



AiF{x,y):-- 



92 



52 



dxf ~^ 8x2 



F{x,y), A2F{x,y):-- 



92 



+ 



dyl dyl 



F{x,y). 



Then, for all x,y,z £ [0,oo)'^, we have [recall (2.1) and (2.3)] 
{z,ViF{x,y)) = {zoy)F{x,y), 
{z,V2F{x,y)) = {zox)F{x,y), 
AiF{x, y) = 8yiy2F{x, y), 
A2F{x, y) = 8xiX2-F(x, y). 



Proof of Proposition 2.1. We use Ito's formula and Lemma 2.2 to 
compute the derivatives of both sides of (2.13) at t = 0: 

^^y[H{{Yt,e),y)]\t=^ 

= {c{9-y),ViF{y,y{l)))F{e,y{2)) 



(2.15) 



and 



+ ^7yiy2AiF(y,y(l))F(0,y(2)) 

H{{y, 0),y)[ci9 -y)o y{l) + ^lymyiWmil)] 



¥.y[H{{y,e),Yt)]\t=o 
= F(e,y(2))f(-cy(l),V2F(y,y(l))) 



d 
di 



(2.16) +^yi(l)y2(l)A2F(y,y(l)) 

+ F{y,y{l)){cy{l),V2F{e,y{2))) 
= H{{y, 9),y)[c{0 - y) oy(l) + 47yiy2yi(l)y2(l)]. 

Since the two derivatives coincide, (2.13) holds (see Corollary 4.4.13 of [6] 
with a = /3 = 0). Equation (2.14) is a direct consequence of (2.13). □ 



Corollary 2.3. Recall Z defined by (2.10). 
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(i) Taking c = 0, Proposition 2.1 implies that Z is self-dual: 

E,[F(Zt,y)] = E,y[F(x,Zt)] for all x,y e [0,oof ,t >0. 

(ii) Letting t — )■ oo ifi (i) (LTid Tccdllijig that /lx\_Zt\ — ^ Qx; 

we get, by 

dominated convergence, the duality relation for the harmonic measure: 



F{z,y)Qx{dz) = / F{x,z)Qy{dz) for all x,y £ [0,oo) . 
Je 

(iii) In particular ( since Qx = Sx for x £ E and due to the symmetry of 
F), for all X £ E and y € [0, oo)^, we have 

[ F{x,z)Qy{dz)=F{x,y) = F{y,x)= [ F{z,x)Qy{dz). 
JE Je 

Corollary 2.4. (i) The family of functions 

J-Q = { [0, oo)2 ^ C : X ^ y), y G [0, oo)2} 

is measure determining for [0, oo)^. 
(ii) The vector space 



(2.17) V:= { }2XmF{;Zm):n£N,Xi,...,XnGC,zu...,ZnGE 




spanned by := {E ^ C:x F(x,z),z £ E} is dense in Ci{E) . In partic- 
ular, T is measure determining for E. 

Proof. Let T>q be the algebra generated by J^o- Clearly, J^o separates 
points of [0, oo)^, contains l = F(-,0) and is closed under multiplication 
and under complex conjugation since F{x, {yi,y2)) = F{x, {y2,yi)). Hence, 
by the Stone-Weierstrass theorem, Vq is dense in the space Ci([0,oo)^) of 
functions [0,oo)^ — )• C that are continuous and have a limit at infinity. As 
J-Q is closed under multiplication, T>q is the vector space spanned by J-q and 
thus J-Q is measure determining on [0, oo)^. 

Let Te = {f\E ■ f £ -^o} ^ and let Ve = {/|_b : / G T^o} denote the alge- 
bra generated by Te- By the above argument, Ve C Ci{E) is dense. Now, 
by Corollary 2.3(iii), an element F{-,y) G Te can be written as the inte- 
gral F{x, y) = J F{x, z)Qy{dz), where the integrand functions are in J-. The 
integral can be approximated (uniformly in x) by finite sums, that is, by 
elements of V . Hence, V is dense in Ve and thus also in Ci{E). □ 

Apparently, Y is ergodic and has a unique invariant distribution with 
a Lebesgue density on (0,00)^. Unlike for the analogous one-dimensional 
equation 



dUt = cib-Ut)dt + y^tdWt, 
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where the invariant distribution is known to be the Gamma distribution 
^20/7,206/7; here, the exphcit form of the density is unknown. It is known (see, 
e.g., [7], Example IV.8.2, page 237) that U hits if and only if 2c6/7 < 1. 
Hence, we may expect that Y = Y^''^'^ hits E only at ((2c02/7) co) x {0}) U 
({0} X (2c0i/7, cxd)). Compare this with the fact that G E' is not hit by the 
infinite 7 process X^'^ (see Theorem 1.7). 

3. Convergence as 7 00: Proofs of Theorems 1.3, 1.4. 

3.1. Construction of the process. Recall the definitions oipt, S and X'^'^ 
in Definition 1.1. In order for the definition to make sense, we still have to 
show, in Proposition 3.2 below, that pt is indeed a Markov kernel and that 
the Chapman-Kolmogorov equation holds. We prepare for Proposition 3.2 
with a lemma. 

Recah the definitions of C, H, D and in (1.8), (1.9) and (1.10). 

Lemma 3.1. (i) D has the Markov property, that is, for x,y £ [0,oo)^ 
and Ac E measurable, we have 

(ii) For / : — 7- C bounded and measurable and r > 0, we have 

/ f{rz)Qx{dz)= I fdQrx- 
Je Je 

(iii) Furthermore, 

I Qx{dz^Qrz+y — Qrx+y 

Je 

Proof, (i) Let F^ denote the filtration generated by the Brownian 
motion B and let F^^ denote the cr-algebra of the Tx past of B [recall (1.7)]. 
Note that F^^dF^. 

For x' G [0, 00)^, denote by P-^' the law of B when started at Bq = —x'. 
Hence, by spatial homogeneity, for x' <x, we have 

P-x' [Br^+y + {x + y)eA] = Qy+^^_^,){A). 
Choosing x' = —B^^ , we infer that 

FBr, [Br^+y + {x + y)eA]= Qy+D^ {A). 
We now apply the strong Markov property of B to obtain 

nDx+y G A I F^] = E[Po[S.,^, + (x + y) G A I F^J, \ F^\ 
= E[PB.jS,^^^ + (rE + 2/)GA]|^f] 



INFINITE RATE MUTUALLY CATALYTIC BRANCHING 



15 



(ii) This follows from the spatial homogeneity of B. 

(iii) Recall that Drx has distribution Qrx- Hence, by (ii) and (i), we get 



Qx{dz)Qrz+y{A) = / Qrx{dz)Qz+y{A) 
JE 



= B[Qy+D.M)]=nDrx+yeA] 

— Qrx+y{A). D 

Proposition 3.2. {St)t>o defined in Definition 1.1 is a Markov semi- 
group. 

Proof. Recall that x i— >• Qx is a continuous map. Hence, for open sets 
A C E, the map x i— t- Qx{A) is lower semicontinuous, by the portmanteau 
theorem (see, e.g., [8], Theorem 13.16), and is hence measurable. Hence, 
X I— > Qx{A) is measurable for all Borel sets Ac E. It remains to check the 
Chapman-Kolmogorov equation for (pt)- By Lemma 3.1 (iii), we infer that 



Pt{x,dy)ps{y,-) = / Qe'"^tx+{l-e-<^t)eidy)Qe-<^''y+il-e-<^'')e 
E JE 

— Qe-':('+»)x+e-<=''(l-e-':*)6»+(l-e-':=)0 

= Qe-<:{t+3)^_^_(^l_^-c(t+3)^g 

= Pt+s{x,-). □ 

3.2. Duality and proof of finite- dimensional distributions convergence (The- 
orem 1.3). In this section, we prove the convergence of the finite-dimensional 
distributions of Y'*''^'^ to those oi X = X'^'^ by means of a duality relation. 
For yi'-'^'^, we have already established the duality, in Proposition 2.1. We 
now come to the duality for X. Recall the definition of Y from (2.11). We 
will need as initial values only y G E x [0,oo)^. Note that, in this case, the 
process Z is constant in time and the process Y is given by the deterministic 
equation 

(3.1) Yt = (e-^*y(l), (1 - e-^*)y(l) + y(2)). 

Hence, Y can be understood as a deterministic Markov process with state 
space E x [0,oof. Recah H from (2.12) and F from (2.2). 

Proposition 3.3. X and Y are dual in the sense that for all initial 
conditions Xq = x £ E, Yq = y £ E x [0, oo)'^ and for all t>0, we have 

(3.2) Bx[H{{Xt,9),y)] = By[H{{x,e),Yt)]- 
In particular, we get 

(3.3) Ex[F{Xt,z)] = F{x,e-'*z)F{e,{l-e-^^)z) for x G [0,oof , z e E, 
and the distribution of Xt is determined by (3.3). 
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Proof. As Y is deterministic, (3.2) and (3.3) are equivalent and so we 
only need to show (3.3). Since z G i?, by Corollary 2.3(iii), the left-hand side 
of (3.3) equals 

/ F{y, z)Q^-ct^+^^_^~ct)Q{dy) = Fie'^^x + (1 - e-'=*)0, z) 

JE 

= F{x,e-^^z)F{9,{l-e-^^)z). 
By Corollary 2.4, equation (3.3) determines the distribution of X^. □ 

We are now ready to prove Theorem 1.3. 

Proof of Theorem 1.3. As both X'^^^ and Y'^''^^^ are Markov pro- 
cesses, it is easy to see that for convergence of finite-dimensional distribu- 
tions, it is enough to show that for any t > 0, x £ E and {x^)^>q in [0, oo)^ 
such that lim^_j.oo x^ ^ x, we have 

(3.4) C,^ [y,^'^'^] [Xf] weakly. 

As shown in the proof of Corollary 2.4(i), Vq is dense in C;([0,oo)^). Hence, 
it is enough to consider F{-,z), z E [0,oo)^, as test functions. Denote by 
Z"' the process defined in (2.10) started at Zq = z. For 7 = 1, we drop the 
superscript, that is, Z := Z^. Denote by the almost sure limit of Zt as 
t — 7- 00 and recall that its distribution is Qz- Note that, due to Brownian 

scaling, {Z^)t>o = {Zyt)t>o- Hence, by Proposition 2.1, we have 

E^^[F{Y,^^''',z)] = E, F{x^,e-''Z^t)F(9,j\e-'''Z^rdr^ 
E,[F(x, e-''Z^)F{e, (1 - e-^*)Zoo)] 
F{x,e-''y)F{e,{l - e-^)y)QMy) 

B,[F{Xt,y)]Qzidy) 



E, 



FiXt,y)Qz{dy) 



= B,[FiXt,z)], 

where the fourth line follows by (3.3) and the last equality follows by Corol- 
lary 2.3(ih). □ 

Remark 3.4. We could also define X'^'^ in Definition 1.1 for initial 
values X G [0, 00)^ (instead of E only). This means that X'^'^ starts life with 
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a jump from x to a random point on E chosen according to Qx and then 
continues with the usual dynamics. Clearly, this process does not have a 
cadlag version (due to the jump at time 0) and its transition semigroup is 
not strongly continuous at 0. Nevertheless, the proof of Theorem 1.3 shows 
that that theorem also holds for this process and hence for Y^''^'^ = X^^ = 
X £ [0,oo)2. 

3.3. Proof of the -convergence (Theorem 1-4) ■ We prepare for the 
proof of Theorem 1.4 with two lemmas. 

Lemma 3.5. Let B = [Bi,B2) he a planar Brownian motion started at 
{Bifi, B2fl) = {u,v) £ [0,oo)^ and let 

T = mi{t>0:Bt^ (0,00)^}. 

Then, for any p € [1,2), we have 

2-p\ttJ 

More generally, one could show for the exit time of a cone with angle 2a 
(here, a = 7r/4) that E[r^/^] < 00 if and only if pa < 7r/2 (see [1], equation 
(3.8)). We give the short proof here in order to be self-contained. 

Proof. By the reflection principle and independence of Bi and B2, we 
get 

P[T>t] = 4Mo,t{0,u)Afo,t{0,v), 

where A/'o,t(a, 6) = (27rt)^"'^/^ e~^^/^* dr is the centred normal distribution 
with variance t. Hence, 

/•oo 

E[rP/2] = / P[T>t'^/P]dt 
Jo 

<f 1 A 0„„*-v.) (i)"\„„,./^. ^ 

Lemma 3.6. For every {u,v) G [0,oo)^, every p G [1,2) and every i = 
1, 2, we have 

I xfQ(„,„)(dx)<|u2-^;2|W2 + ^!{!M!^<^. 
Je cos(p7r/4) 
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Proof. This can be verified by means of an explicit computation using 
the density formula of Q(^u,v) '^^ (2.5). □ 

Note that finiteness of the expression on the left-hand side in Lemma 
3.6 (which is what we need in the proof of Theorem 1.4) could also be 
inferred without computations by the Burkholder-Davis-Gundy inequality 
and Lemma 3.5. 

Proof of Theorem 1.4. (i) By Lemma 3.6, we have 

E[{X^fr]= I yfQ,-c.,+(i_,-c.),(dy)<oo. 
JE 

Fix t > and define 

Ml := e-^'x, + (1 - e-'')9i + e<^'-'^ \l lY^fY^f dW,,^. 

Let {M\) = (Ml) denote the square variation process of both M\ and Mg. 
Note that M*j = Y^f'^ > and that M* is a martingale and thus 

(3.5) Ml = I Ml] > for all s G [0, t]. 

Now, (M*) s>o is an isotropic diffusion in and is hence a time-transformed 
planar Brownian motion. That is, there exists a planar Brownian motion B 
(with respect to some right-continuous complete filtration J^) started at 
Bq = e~'^^x + (1 — e~'^^)6 such that each {M\ ,)s is an F stopping time and 
such that -B^jvft = M* for all s > 0. 

Define the J- stopping times 

T ■.= mi{s> Q:Bsi (Q, oof ] and tq := inf{s > : ^ [0, oo)^}. 

Clearly, we have r = tq almost surely and, hence, by (3.5), 

{Ml)t<TQ = T a.s. 

Using the Burkholder-Davis-Gundy inequality for the martingale (i?i,s)s>o 
yields (see Lemma 3.5) 

EfsupS^? 1 < 2'P-\Bl^ + (4p)PE[rP/2]) < oo. 

Hence, (|-Bj,rAs|^)s>o is uniformly integrable and we can apply the optional 
sampling theorem to obtain 

(ii) By Theorem 1.3 and the Skorohod embedding theorem, we may con- 
struct all processes on one probability space such that Y^"'^' — )• X^' almost 
surely as 7 — )• oo. By part (i), the pth moments of Y^^^' , 7 > 0, are uniformly 
integrable and so we have the desired L^-convergence. □ 
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4. The strong construction (proofs of Theorems 1.6 and 1.7). Recall the 
definitions of C, H and D in (1.8) and (1.10). 

Lemma 4.1. The map x^Dx is cddldg. 

Proof. This follows from continuity of B and the definition of r^. □ 

Proof of Theorem 1.6. From Lemmas 3.1 and 4.1, we infer that X'^'^ 
has the Markov property and cadlag paths. It remains to show (1.11). 
By Lemma 3.1, for x,z £ E, Ac E measurable and < s < t, we have 

c 6 

(with Px denoting the probability law of X^' , as defined in Theorem 1.6) 
Ps,t{z,A) = Px[xf G A \ Xf = z] 

= P[C(0, t)Dx+E^o,t) e A I I),+s(o,s) = ^^(0, syh] 

= Qc{0,s)-iz+Eis,t)(.C{^,t)~^A.) 

= Qc{s,t)z+CiO,t)E{s,t)iA)- □ 

Proof of Theorem 1.7. U c9 = 0, then X"^'^ is the deterministic pro- 
cess X^'^ = c'^^Xq^ and hence is polar. 

Now, assume that c6 ^ 0. Le Gall and Meyre [12] show that almost surely, 
for all z £ (0,oo)^, the planar Brownian motion B does not leave the cone 
[—z,oo) first at —z. More formally, consider the event 

A := {Br, / -z for ah z £ (0,00)^}. 

Theorem 1 of [12] then implies that P[^] = 1 (in fact, they show that no 
rectangular cone is first left at its vertex, not only northeast cones [z,oo)). 
Now, by (1.12), we have 

{X^'^ / for aU t > 0} = {D^+re / for all r > 0} 

= {Br^_^^g ^x + rO for all r > 0} D A. 
This shows the claim of Theorem 1.7. □ 

5. The infinitesimal dynamics of X^'^. In this section, we give a de- 
scription and construction of the infinite rate mutually catalytic branching 
process X in terms of its infinitesimal characteristics. To this end, we will 
define a linear operator Q'^'^ that: 

(i) defines the contraction semigroup of X in the sense of the Hille- 
Yoshida theorem (Theorem 5.3); 

(ii) defines a well-posed martingale problem whose unique solution is X 
(Theorem 5.4). 



20 A. KLENKE AND L. MYTNIK 

5.1. Results. Recall, from Definition 1.1, that the linear operator St on 
Ci{E) is defined by 

Stf{x):= I f{y)pt{x,dy)= / f{y)Qe-'=tx+{i-e-'=t)eix,dy). 

J E J E 

In order to define the generator of 5 = {St)t>o-, we will need to study (for 
suitable functions /) the limit 

(5.1) limt-\Stf{x)-f{x)) = hme-'(^l f dQ,+,,^e-.) - fi^)^- 

In the sequel, we will use the shorthand notation 

d d 
dif{u,v):=—f{u,v) and d2f{u,v):=—f{u,v). 

In order to define what we mean by a suitable function, we introduce the 
subspace Cf{E) C Ci{E). 

Definition 5.1. Let Cf{E) c Ci{E) be the subspace of such functions 
f^Ci{E): 

(i) whose partial derivatives dif and ^2/ exist on (0,oo) x {0} and 
{0} X (0,00), respectively, are continuous, can be continuously extended to 
{0} X [0, 00) and fulfill 

(5.2) lim nai/(n,0)= lim vd2f{Q,v) = {); 



v—^oo 

(ii) whose partial second derivatives dff and (9|/ exist on (0, 00) x {0} 
and {0} X (0,oo), respectively, and are such that 

(5.3) 11/112,00 := sup r(|a2/(r,0)| + |a|/(0,r)|)<oo. 

re{0,oo) 

Note that, for / e Cf{E), we have 

(5.4) ||/||i,oo:= sup (|ai/(r,0)| + 1^2/(0,01) < 00. 

re[0,oo) 

In order to get an explicit formula for the limit in (5.1), we define the 
vague limits (for u, v > 0) 

V{}d,v)-=v-J^^£~^Q(e,v) and z/(„^o) :=^'-lime~"^QKe)• 
^(u,o) can be thought of as the "Levy measure" of the next jump when the 
actual position is (u, 0) and the drift is in direction of (0,1). In order to 
formalize this, for the drift in direction (0,1), we define the linear operator 
02 on Cf{E) by ^2/(2:) = d2f{x) if xi = and 

Q2f{x) = [ [f{y) - fix) - (yi - xi) dif{x)]u^{dy) if xi > 0. 
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For the drift in direction (1,0), we define Gi similarly. Note that Vx is not 
a finite measure and that the integral of yi — it with respect to i^(n,o) is well 
defined only as a Cauchy principal value and, as such, equals zero. Hence, 
this term in the integral is needed in order for the integral to be well defined 
in the usual sense. In Lemma 5.5 below, we will show that Qif and ^2/ are, 
in fact, well defined and are in Ci{E). 

Due to spatial homogeneity of planar Brownian motion, we have a scaling 
relation that helps to get rid of the many different I'x in the definition of Qi 
and Q2: 

f{x)i'(u,o){dx) = ^ j^f{ux)v(ifl){dx). 
Furthermore, letting /''"((xi, X2)) := f{{x2,xi)), by symmetry, we have 

f{x)vfs),v){dx)= I P{x)u(^y^o){dx) = - I f\vx)u(^ifi){dx). 
Je V Je 

Hence, we can express Gi and G2 in terms of 
(5.5) i/:= 1/(^0). 

Using the explicit form of the density of Q[i^£) in (2.5) and letting e — t- 0, 
we get that the a-finite measure u on E has a one-dimensional Lebesgue 
density given by 

(A u 



vr (1 — ?i)^(l + li)^ 
4 V 

dv, if « = 0. 



(5.6) u{d{u,v)) = 

Qi and Q2 can now be written as 

{d2fix), ifxi = 0, 

- / ifixw) - fix) 
xi Je 
-xi{yi -l)di f{x)]iy{dy), if xi > 0, 

and 

(5.8) gif = {G2p)^. 

Finally, we define the operator ^^'^ on Cf{E) with domain 'D{Q'^'^) = Cf{E) 
that determines the infinitesimal characteristics of the process X = X'^'^: 

2 

(5.9) g'''f{x) = ^c{ei-x,)gif{x). 

i=l 
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Lemma 5.2. The operator G^'^ is well defined. That is, for f G Cf{E), 
the expressions in (5.9) and (5.7) are well defined and we have G'^'^f G 
Ci{E). 

This lemma will be proven in Section 5.2. 

Theorem 5.3 {X^'^ via its generator), (i) For every f G Cf{E), we 
have, pointwise, for all x & E, that 

(5.10) g'''f{x) = \\uie~^( f fdQ^+,,^e-.)-fix))=\im^^i^^^f^. 

>^i'^ \Je J '-l-O * 

(ii) The operator Q'^'^ on Ci{E) is closahle and its closure generates the 
contraction semigroup S of the process X'^'^ . 

The theorem will be proven in Section 5.2 using the classical Hille-Yoshida 
theorem. 

A different, and more modern, approach to constructing Markov processes 
from their infinitesimal dynamics is the martingale problem technique due 
to Stroock and Varadhan. 

Recall from (2.17) that V C Cf{E) is the vector space spanned by {F{-,z), 
z G E^. Define the linear operator Q'^'^ on V by (5.9) and (5.7). By Theorem 
5.3(i), we obtain for z £ E [using Corollary 2.3(iii) in the second line and 
Lemma 2.2 in the last line] that 

g'''F{-,z){x) = hme-i ( [ F{y, z) dQ,+,,^e-x){dy) - F{x, z) 

= \im.£~'^{F{x + ec{e - x),z) - F{x,z)) 

(5.11) 

= {c{d-x),ViF{x,z)) 
= F{x,y)[c{9 — x)o z]. 

Hence, (5.11) is enough to define G'^'^ on V and we do not really need the 
measure v from (5.7) here. 

A solution of the {G'^'^, V) martingale problem is an valued measurable 
stochastic process X such that 

Mt:=F{Xt,z)- [ {c{9-Xs)oz)F{Xs,z)ds 
Jo 

is a (C-valued) martingale. A martingale problem is said to be well posed 
if, for every probability measure fj, on E, there exists a solution X with 
£[Xo] = /Li (existence) and any two solutions have the same finite-dimensional 
distributions (uniqueness). In this case, X is a Markov process (see [6], 
Theorem 4.4.2(a)). 
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Theorem 5.4 (Martingale problem characterization of X'^' ). The mar- 
tingale problem {G'^'^,V) is well posed and its unique solution is X'^'^ . 



This theorem will be proven in Section 5.3. 



5.2. The Hille-Yoshida approach (proof of Theorem 5.3). Lemma 5.2 
and part (i) of Theorem 5.3 are direct consequences of the following two 
lemmas. 

Let ei = (1,0) and 62 = (0,1). 

Lemma 5.5. For f G Cf{E), x £ E and i = 1,2, the expression Gif{x) 
from (5.7) and (5.8) is well defined and we have 

(5.12) hme-i ( I fdQ,+ee, - fix)) = Gifix). 



Lemma 5.6. For f G Cf{E), we have G"'^ f G Ci{E). 



Proof of Lemma 5.5. For x = (0,0), since Qsa = ^eci, this is the very 
definition of Gi- For n 7^ (0,0), by linear scaling and symmetry, it is enough 
to consider the case x = (1,0). If i = 1, then the left-hand side of (5.12) 
equals 

lime-i(/(l + e,0) - /(1,0)) = 5i/(l,0) = (gi/)(l,0). 

£4-0 

Now, consider i = 2. It is a simple exercise to compute that for every e > 0, 
4 r°° r(r-l) 2 



vr 



/■°° r(r- 1) , 2 . 

I , o , / 9 , — 9 — TTo dr = -e arctan(e) 
Jq 4e2 + (r2 j^gi -ly TT 



4 



ds. 



vr 7o 4e2 + (^2 _ ^2 + 1)2 
Hence, if we let g{y) := (yi — 1) 5i/(l, 0), then, for every e > 0, 

(5(y)-5(l,0))Q(i,e)(fiy) = 0. 

Hence, we can replace f hy f — g. Now, f — g is twice differentiable, has at 
most linear growth and di{f — g){l,0) = 0. Hence, 



|(/-g)(n,0)- /(l 

u>0,u^l (u - 1) 



sup ;^^2 < 
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This allows us to use dominated convergence in the following computation 
to obtain 

lime-i(^^/ciQ(i,o)+ee.-/(l,0)) 

= lime-i / [/(x)-/(l,0)-(rEi-l)5i/(l,0)]Q(i,,)(dx) 



4e2 + („2_e2 + 1)2 

liml^; „ . "'if'"''"'"''/':-"'' „ 

'2e^{u + l){u-l) + {u + l)^{u-l)^ 

4 r-- v[{f-g)(0,v)-f(l,0)] 

4 r ^[(/-9)(^,o)-/(i,o)]^^ 



'4 




vr 


/ 




4 


+ 






TT 


'4 




vr 





vr7o (^2-1)2 

4 r-.[(/-,)(0,.) /(1,0)]^^ 



[/(y)-/(l,0)-(2/i-l)9i/(l,0)]i.(dy) 

= ^2/(1,0). □ 

Proof of Lemma 5.6. We have to show that for any / G Cf{E), 
Q'^'^ f{x) is continuous in x G and has a limit at oo. By (5.9), it is enough 
to derive these properties for Gi{x) := {9i — Xi)Qif{x), i = 1,2. We will give 
the proof only for the case i = 2 since the case i = 1 is analogous. 

For xi = 0, we have 

(5.13) G2{x) = G2(0, X2) = {62 - X2) 52/(0, X2). 

This expression is clearly continuous in X2 G [0,oo) and, by (5.2), we have 

(5.14) lim G2(x) = 0. 

Now, consider the case xi > 0. Hence, by (5.7), 

G2{x) = j g{x,y)u{dy), 

where 

g{x,y) := —[f{xiy) - f{x)-xi{yi - l)dif{x)]. 

Xl 
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Since / G Cf{E), for all y ^ E, we have: 

(i) X ^ g{x,y) is continuous on (0,oo) x {0}; 

(ii) limj;^^oogix,y) = 0; 

(iii) g{x,y) = 6*2 ^2/(0, 0)y2- 

In order to find an integrable dominating function for g, define h:E^ [0, 00) 
by [recall (5.3) and (5.4)] 



h{y) :-- 



'2IU ||2,oo(?/i - l)^ ifyiG(i,|), 
2^211/111,00(2/1 + 2/2 + 1), otherwise. 



Note that the density of i^idy) decays like 1/(2/1 + 2/2)'^ as y — )• 00. Further- 
more, 

( -|a2 '^(^2/) ^ 4 yi 

' dyi 7r(l + yi)2 

is bounded on (1/2,3/2) x {0}. Hence, we have f hdv < 00. 
For all y € £^ and xi > 0, we have 

\g{x,y)\ < -ilfixiy) - /(0,0)| + |/(x) - /(0,0)| + xi(yi + l)|5i/(x)|) 

Xi 

< 2^2(2/1 + 2/2 + 1)11/111,00. 

Furthermore, recalling (5.3), for yi G (1/2,3/2), by Taylor's formula, we get 
that 

\g{x, y)\ = — |/((2/i - l)x + x) - f{x) - xi{yi - 1) dif{x)\ 

Xl 

<^(yi-l)2 sup xi\dlfiu,0)\ 

^ u>xi/2 
< ^211/112,00(2/1 -1)'. 

Hence, in fact, \g{x,y)\ < h{y) for all y £ E, x £ (0,oo) x {0} and the domi- 
nated convergence theorem yields that G2 shares the properties (i) and (ii) 
of g{x,-) and that 



limG2(x) = 02^2/(0,0) / y2Z^(d2/) = ^2^2/(0,0) = G2(0,0). 

xi^O J 

Combining this with (5.13) and (5.14), we have G2 G Ci{E). □ 

In order to show part (ii) of Theorem 5.3, we will apply the Hille-Yoshida 
theorem for generators of contraction semigroups. Recall, from Corollary 
2.4, that V is dense in Ci{E). Also, by Lemma 2.2, one can easily check that 

VdCfiE). 
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For each z G E, define the map Uy : [0, oo) — Ci{E) by Uy{t) := StF{-,y). 

By [6], Proposition 1.3.4, the operator Q'^'^ on Ci{E) with domain 'D{Q'^'^) = 
Cf{E) is closable and its closure generates (uniquely) the semigroup (5i)t>o 
on Ci{E) if the following conditions are all fulfilled: 

(a) Q'^'^ is dissipative; 

(b) Uy{t)£V{g<'^^) for ah t>0; 

(c) the map (0, oo) -^Ci{E), t^Q'^'^Uy{t) is continuous; 

(d) for ah t > 0, 

(5.15) Uy{t)-Uy{Q)= [ g'''%y{s)ds. 



Hence, in order to prove part (ii) of Theorem 5.3, it remains to check (a)-(d). 

(a) Let / € Cf{E) and assume that / assumes its maximum at x S 
E \J {oo}. Since Stf{x) < f{x) for all t > 0, equation (5.10) implies that 
Q'^'^ f{x) < 0. Hence, Q'^'^ fulfills the positive maximum principle and is thus 
dissipative (see, e.g., [6], Lemma 4.2.1). 



(b) By Proposition 3.3, for any y £ E, x £ E and t > 0, we have 
(5.16) uy{t){x) = StF{;y){x) = F{x,e~''y)F{9, (1 - e"^*)?/). 

As F(-,e-^*y) is in Cf{E), so is StF{;y). 

(c) By (5.10), we have 

g^'^uymx) = \ime~\St+e - St)F{.,y){x) = ^{uy{t){x)). 
£\Q at 

Using (5.16) and Lemma 2.2, for every x £ E, we get 



implies (5.15). 

5.3. The martingale problem (proof of Theorem 5.4)- Before we prove 
this theorem, we derive a duality relation for processes satisfying the mar- 
tingale problem {G'^'^,V). Recall the definition of Y from (2.11). 
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Lemma 5.7. Let fi be a probability measure on E. Let X be any solution 
of the martingale problem {Q'^'^ ,V) with C[Xq\ = /i. Then X and Y are dual, 
in the sense that for any y £ E x [0,oo)^, we have 

E^[H{{Xt,e),y)]= [ By[H{{x,6),Yt)]fi{dx) forallt>0. 

JE 

Proof. As X is a solution of the martingale problem, we have that 

H{{Xt,9),y) - f H{{Xs,9),y)[c{9 - Xs) o y{l)] ds 
Jo 

= F{9, y{2)) (^F{Xt,y{l)) " F{Xs,mM9 - X^) o y{l)] ds^ 

is a martingale. On the other hand, by (2.16) [since y{l) G E, one term 
vanishes] , 

jEy[H{{x,9),Yt)]\t=o = (-cy(l), V2F(x,y(l)))F(e,y(2)) 

(5.17) + (cy(l), V2F(0,y(2)))F(x,y(l)) 

= H{{9,x),yM9-x)oy{l)]. 
Since Y is deterministic, we get that 

H{{x,9),Yt) - f H{{9,x),%)c{9 - x)oYs{l)ds = H{{x,9),y) 
Jo 

is the trivial martingale. By [6], Corollary 4.4.13, this implies that 

B^[H{{Xt,9),y)]= J Ey[H{{x,9),Yt)Mdx) 
and we are done. □ 

Proof of Theorem 5.4. By Theorem 5.3(ii) and (5.11), and since 
V C Ci{E), by definition of X^'^, the process X'^'^ is, in fact, a solution of 
the martingale problem (G'^'^jV). 

Now, assume that X and X' are two solutions with i2[Xo] = ^[^q] ~ ^• 
By Lemma 5.7, we get 

B^[F{Xt, y)] = E^[F{Xl, y)] for ah t > and yeE. 

By Corollary 2.4, {F(-,y),y S E} is measure determining on E. Hence, 
Cf,[Xt] =Cf,[Xi] for ah i > 0. By [6], Theorem 4.4.2, this implies that the 
finite-dimensional distributions of X and X' coincide. □ 
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